HIGH ORDER MULTI-SCALE WALL-LAWS, 
PART I : THE PERIODIC CASE 
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Abstract. In this work we present new wall-laws boundary conditions including microscopic os- 
cillations. Wc consider a newtonian flow in domains with periodic rough boundaries that we simplify 
considering a Laplace operator with periodic inflow and outflow boundary conditions. Following 
the previous approaches, see [A. Mikelic, W. Jager, J. Diff. Eqs, 170, 96-122, (2001) ] and [Y. 
Achdou, O. Pironneau, F. Valentin, J. Comput. Phys, 147, 1, 187-218, (1998)], we construct high 
order boundary layer approximations and rigorously justify their rates of convergence with respect 
to e (the roughness' thickness). We establish mathematically a poor convergence rate for averaged 
second-order wall-laws as it was illustrated numerically for instance in [Y. Achdou, O. Pironneau, F. 
Valentin, J. Comput. Phys, 147, 1, 187—218, (1998)]. In comparison, we establish exponential error 
estimates in the case of explicit multi-scale ansatz. This motivates our study to derive implicit first 

- 

order multi-scale wall-laws and to show that its rate of convergence is at least of order e 2 . We provide 
a numerical assessment of the claims as well as a counter-example that evidences the impossibility 
of an averaged second order wall-law. Our paper may be seen as the first stone to derive efficient 
high order wall-laws boundary conditions. 
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layers, finite element methods, error estimates. 
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1. Introduction. The main goal of wall-laws is to remove the stiff part from 
boundary layers, replacing the classical no-slip boundary condition by a more sophis- 
ticated relation between the variables and their derivatives. They are extensively 
used in numerical simulations to eliminate regions of strong gradients or regions of 
complex geometry (rough boundaries) from the domain of computation. Depending 
on the field of applications, (porous media, fluid mechanics, heat transfer, electro- 
magnetism), wall-laws may be called Beavers- Joseph, Saffman- Joseph, Navier, 
Fourier, Leontovitch type laws. 

High order effective macroscopic boundary conditions may also be proposed if we 
choose a higher degree ansatz, see [7] for applications in microfluidic. In a similar 
perspective but in the context of fluid mechanics, numerical simulations have shown 
that second order macroscopic wall-laws provide the same order of approximation as 
the first order approximation. Recently a generalized wall-law formulation has been 
obtained for curved rough boundaries |17[ I19| and for random roughness [4] . Note 
that such generalizations are important from a practical point of view when dealing 
with e.g. coastal effects in geophysical flows. From a mathematical point of view, 
wall-laws arc also interesting. In the proof of convergence to the Eulcr equations, 
the 2D Navier-Stokes system is complemented with wall-laws of the Navier type [5]. 
Recently several papers analyze in various settings the properties of such boundary 
conditions, see [12], [16], [U], [S], [13]. 

In this paper, we focus on fluid flows. Starting from the Stokes system, we sim- 
plify the problem by studying the axial velocity through the resolution of a specific 
Poisson problem with periodic inlet and outlet boundary conditions. Our scope is 
to justify mathematically higher order macroscopic wall-laws and to explain why in 
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their averaged form they do not provide better results than the first order laws. We 
shall explain how to get better estimates including some coefficients depending on the 
microscopic variables: this leads to new oscillating wall-laws. 

The basic scheme to establish standard averaged wall-laws is the following (see fig. 
II. First we use an ansatz for the velocity and the pressure which will give, after an 
adequate extension, a main order term completed with some boundary layer correc- 
tors defined on the whole rough domain (fig. 1.1, step I). This is possible due to the 
boundary layer theory that can be seen as a particular case of a general homogeniza- 
tion process. In a second time, a specific average is performed on this approximation 
and a new boundary condition of mixed type is recovered on a smooth fictitious inter- 
face strictly contained in the domain (fig. 1.1, step II). As one sees on the figure the 
only difference between Achdou's and Jager's approaches is situated in the boundary 
layer's construction. It is an easy task to show that they are in fact a specific lift one 
of the other. 
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figure 1.1: The standard approach: from the exact solution to an averaged wall-law 



The main result in our paper is the derivation of a high order boundary layer approxi- 
mation that satisfies the homogeneous Dirichlct boundary condition on the rough wall 
and that leads to new wall-laws with microscopic effects see fig. 11.21 The ansatz is 
expanded up to the second order in e and an exponential convergence in the interior 
domain is obtained using it, (fig. 11.21 step I'). Despite this great rate of convergence, 
the corresponding second order averaged wall-law behaves badly and does not preserve 
the nice convergence properties of full boundary layer approximations. The estimates 
show that this is due to the great influence of microscopic oscillations. We then de- 
rive new wall-laws that do converge exponentially on the smooth domain. They have 
the form of explicit non-homogeneous Dirichlet boundary conditions and they depend 
on the zcroth order Poiseuille flow as well as on the microscopic oscillations on the 
fictitious interface (fig. 11.21 step IF). 

At this stage, we go one step further and derive an implicit multi-scale first order 
wall-law. We obtain a Saffman-Joseph's like law that now contains a coefficient 
that includes the microscopic oscillations. We rigorously derive a rate of convergence 
in e5, thanks to the steps introduced in the previous sections (fig. If .21 step II"). 
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figure 1.2: The new approach: from the exact solution to multi-scale wall-laws 



We underline that this work is a necessary building block when studying wall-laws 
for the stationary Navier-Stokes equations: asymptotic expansion of the quadratic 
non-linearity transfers a cascade of contributions to the microscopic cell problems, as 
already noticed in [2]. The first order cell problem is homogeneous and the second- 
order cell problem involves the non-linearity of the first order approximation. Until 
now, every averaged wall-law was only first order accurate and thus wall-laws were 
not able to display second order effects of non-linearities. 

In a wider context that does not concern only fluid flows, the main concept this 
work emphasizes is the following: we have shown that it is possible to replace a ge- 
ometrical roughness and "smooth" boundary conditions (in the sense unperturbed, 
as for instance homogeneous Dirichlet ones) by a smooth domain but with a multi- 
scale perturbed boundary conditions, (see fig. 11.31 below). Depending on the kind of 
boundary perturbation, we get different orders of precision in this process. For com- 
plex multi-scale 3D problems, we still expect some numerical gain when performing 
this switch, especially if one uses some increased multi-scale finite element bases (see 
[9] and references therein). 
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figure 1.3: One of the main points of this article: switching perturbations from geom- 
etry to boundary data. (BC stands for boundary conditions) 



To show the practical importance of the results above, in Section [6l we perform 
numerical tests on a 2D case. For various values of e, we first compute the rough 
solution it^ on the whole domain fl'^ , then we compute the wall-law solutions defined 
only on the interior smooth domain . We perform these tests in the periodic case. 
We recover exactly theoretical claims: numerical error estimates confirm that averaged 
wall-laws do not differ at first and second orders. We prove that our new implicit multi- 
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scale wall-law provides better results than classical averaged laws. However, the fully 
explicit approximations still show higher order convergence rates with respect to e. 

2. The simplified problem: from Navier-Stokes to Laplace equation. 

In this work, il'^ denotes the rough domain in M? depicted in fig. 12. 1[ denotes 
the smooth one, is the rough boundary and F*' (resp. F^) the lower (resp. upper) 
smooth one (see fig. 12. 

Hypotheses 2.1. The rough boundary is described as a periodic repetition at the 
microscopic scale of a single boundary cell . The latter can be parameterized as the 
graph of a Lipshitz function f : [0, 2t:[—>- [—1 : 0[ such that 

P" = {ye [0,27r] X [-1 : 0[/y2 = fivi)} (2.1) 

Moreover we suppose that f is negative definite, i.e. there exists a positive constant 
5 such that f{yi) < 6 for all yi E [0, 2tt]. We assume that the ratio between L (the 
width of ilP) and 27re (the width of the periodic cell) is always an integer called N. We 
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figure 2.1: Rough, smooth and cell domains 



consider a simplified setting that avoids the theoretical difficulties and the non-linear 
complications of the full Navier-Stokes equations. Starting from the Stokes system, 
we consider a Poisson problem for the axial component of the velocity. The pressure 
gradient is assumed to reduce to a constant right hand side C. We consider only 
periodic inflow and outflow boundary conditions. The simplified formulation reads : 
find u*^ such that 

- Au" = C, for X e n', 

u'=0, xeruT\ (2.2) 
u*^ is Xl periodic. 

We underline that the results below can be directly extended to rough domains with 
smooth holes and to the Stokes system. 

In what follows, functions that do depend on y ~ x/e should be indexed by an e (e.g. 
Ue ~Uc{x,x/t)). 

3. The full boundary layers correctors. 

3.1. A zeroth order approximation. When e = 0, the rough domain Q,^ 
reduces to which is smooth. The solution of system (|2.2p in this limit is known 
and explicit: it is the so-called Poiseuille profile : 

u''{x)^^{l^X2)x2, Va-er!°, 
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the latter term should be our zeroth order approximation when performing an asymp- 
totic expansion w.r.t. e for e > 0. The determining step is then how to extend this 
zeroth order approximation so that it is defined on the whole domain Q"^. A possible 
choice is to use the Taylor expansion of near X2 = 0, it leads to define the zeroth 
order expansion as a C^{fl'^) function that reads 

(u°{x), if.Tef]° 
y. dx2 

Remark that this particular choice does not satisfy the homogeneous Dirichlet bound- 
ary condition on F"^. Next we estimate the zeroth order error wrt the exact solution. 
Proposition 1. If is a open connected piecewise smooth domain, the solution 
exists in H^(fl'^) and is unique. Moreover we have 



where the constants Ci and C2 are independent on e 

Proof. It is based on standard a priori estimates and a duality argument. The 
existence and uniqueness of u*^ are standard and left to the reader. We focus only on 
the error estimates. Namely, r° := — u° satisfies 



r 







Ar" ^ Cx[ir\no], in 0% 
0, on r\ 

r" ^~j—{xi,0)x2 on T° , 

OX2 

is xi — periodic on Fin U Fout 

There, one remarks that a part of the error comes from the source term localized in 
fJ*^ \ ri°, and another part comes from the non homogeneous boundary term on F*^. 
We set the lift 



ou 
dx2 



X2X[Q.'\iV'], and z := r° 



then the weak formulation reads : 

(V2, Vt;)n. = {C,v)n>\no - (C,zOro, v e H^^i^'), 

where the last term in the rhs comes when applying the Laplace operator A on s. 
Thanks to Poincarc-like estimates we have the following properties of the norm 
and the semi-norm on fi*^ \ fi" 




n'\n° 




< C5(e + Ve)y J^^^^^^^^ |Vif ) < 2c5V^||v||si(o.). 

This leads to the H^{il'^) estimate. For the norm, we use the concept of a very 
weak solution Namely, one solves the dual problem: for a given ip e L^(f7"), ip 
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being xi — periodic on Fin U Tout find v e H'^{n^) such that 

v = o, Va;er°ur\ 

w is xi — periodic on Fin U Tout- 
Considering the L^(ri°) scalar product, and using the Green formula 

\9n /g^o \dn Jq^^o 



(3.1) 



where the brackets refer to the dual product in iJ 2 (r°), and the rest of products are 
in 1? ^ cither on F*^ or on 17°. Then, one computes 



( '^'^ II II II nil 

Vail'" Jro ^^6||^IL.(no)|k IL.(ro)- 



The last estimate is obtained thanks to a linear dependence of the normal derivative 
of the trace of v on the data <~p^ HH]. Thanks to Poincare estimates, one writes 

which ends the proof by taking the sup over all in □ 

3.2. A first order correction. The zeroth order correction contains two dis- 
tinct sources of errors : a part is due to the order of the extension in fJ"^ \ and 
another part comes from a non homogeneous rest on . In what follows we show 
that a first order extension can be corrected by series of terms that makes the full 
boundary layer approximation vanish on . 

The micrscopic cell problem : In order to correct on V^, one starts by solving a 
microscopic cell problem that reads : find (3 s.t. 

- A/3 = 0, in Z+ U P, 
/? ^ -y2, on P", (3.2) 
/? is Hi — periodic . 

We define the microscopic average along the fictitious interface F : 

As Z+ U P is unbounded in the 1/2 direction, we define 

D^'^ = {ve LlciZ^ U P)/Dv e L^{Z+ U P)^ visyi- periodic }, 
then one has the following result : 

Theorem 3.1. Under hypotheses \2.1l there exists (3, a unique solution of p.2p 
belonging to D^'^ . Moreover, there exists a unique periodic solution rj E i72(F), of 
the following problem 

< S'77, ^ >=< 1,/i >, V/iei7^(F), 
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where <, > is the _ff"5(r) — H^(r) duality bracket, and S the inverse of the Steklov- 
Poincare operator (see appendix \A.l\) . One has the following correspondance between 
(3 and the interface solution rj : 

where Hz+rj (resp. Hprj) is the yi-periodic harmonic extension of rj on (resp. 
P). The solution in can be written explicitly as a series of Fourier coefficients of 
rj and reads : 

r2TT 



oo 

k — — QO 



r^{yi)e-'''y'dy,. 



In the macroscopic domain fi" this leads to 



< K^/^\\n\\ 1 . 



The proof is given in the appendix for sake of conciseness. The corresponding macro- 
scopic fuU boundary layer corrector should contain at this stage 



where we subtract /? in order to cancel /?'s errors on . 

In order to cancel the contribution of the constant (3 near the rough boundary but 
keep its benefit close to F^, one solves the "counter-flow" problem: find d s.t. 



-Ad = 0, mn°, 

1 on r°,(i on r\ 
d is xi — periodic on Fin U T^utj 



(3.3) 



the solution is explicit and reads = (1 — X2)- Moreover, it can be extended to the 
whole domain O'^. The complete first order approximation now reads : 

:= ^0 + o)(/3 - /3) + e^(xi,0)/3(l - x^), Vx £ 

0x2 0x2 

= < + e^(xi,0)(/3-/?X2), 
0x2 

the first index of u^^^ corresponds to the extension order of vP in Qf" \ £7", while the 
second index is the order of the error on F^. Indeed, if we consider the trace of u\''^ 
on F*^, we have a second order error 



,1.2 I 



2 ( (^^l-o\ ^2 



dx: 



dul- 



dx2 



Again, this error is linear and should be corrected by the micro boundary layer (3. 
A similar macroscopic boundary layer correction process should be performed at any 
order leading to 

-^"1 ( p^^^ 



e-^(xi,0) 
0x2 



PX2]+€I3\I3 



13X2 



-e^f^3[^^-T3x2)+. 



1 + e/? dx: 



e 

I3X2 



(3.4) 



This approximation satisfies a homogeneous Dirichlet boundary condition on F*^, and 
solves 



~ Au]^^ ^ Cx[n^\no], in 
ul^°^ = 0, on ^^ 



is xi — periodic on Fin U Fout- 

If we consider the corresponding approximation error, we obtain 

Proposition 2. Under hypotheses \2.1l the error of the first order approximation 

satisfies 

- ul-°°\\Hiin') ^ cse, II"' - w.''°°|L2(oo) < cge^ 
where the constants csjCg are independent on e. 

The proof follows the same lines as in proposition[l]except that the significant source of 
errors is the rhs of the first equation in (|3.5[) , while an exponentially small microscopic 
perturbation lies on F^, on the contrary there are no errors on F*^, because = 
on it. 

3.3. Second order approximation. Instead of extending only linearly the 
Poiseuille profile it is obvious that a quadratic term is missing to complete the ap- 
proximation. In the following denotes the second order extension of iP in 51^ \ 51°. 



(J 



^^_^-(a;i,0)a;2 + ^^(a;i,0)^, x e X[n^\n°] ^ 

The second order error on F*^ is corrected thanks to a new cell problem : find 7 e _D^'^ 
solving 

A7 = 0, in Z+ U P, 
7 = -vl. on po, (3.6) 
7 periodic in y\. 

The proof of the following proposition is left in the appendix IA.2I 
Proposition 3. Under hypotheses \2.1[ there exists a unique solution 7 0/ p.6p 
in D^'^^Z'^ U P). Moreover it admits a power series of Fourier modes in and 
7 e [-1,0] if Pa [0,27r] X [-1,0]. 

The horizontal average is denoted 7. The same multi-scale process leads to write the 
full boundary layer approximation as 

Again a third error remains on F'^ and it is linear wrt to 2/2, thus it should be corrected 
thanks to the series of first order cell problems as in (|3.4p . We set u'^'°° to be the second 
order approximation that satisfies a homogeneous Dirichlet boundary condition on F*^, 



it reads 



l + epdx2 V V e / 



+ 



2 dx: 



'-{xuO) 



Our approximation satisfies the following boundary value problem 
f - Auf°° = C, in 
= 0, on r% 

^ is xi — periodic on Fin U Tout, 
where g is the contribution of the microscopic correctors on and reads ; 

1 r^n n r 

67 



(3.7) 



2 a^(^^'O) 



7(^,1) -7 



Remark that the only error remains on and as the proposition below claims, it is 
exponentially small wrt e. 

Proposition 4. Under hiivotheses \2. 1\ the error of the first second order approxima- 
tion satisfies 



2,oo|| ^ /- -- 



where the constants cg, cy are independent on e. The proof is identical to the one of 
proposition [T] except that the only source of errors is the contribution of function g^, 
there arc nor errors on F"^, neither source terms inside fi"^. 

4. Averaged wall-lav^rs. 

4.1. The averaged v^rall-lav^fs: a nev^r derivation process. At this stage, we 
rewrite our first and second order approximations separating slow and fast variables 



,l,oo _ „,0 



1 + e/3 aa;2 1 + e/3 ^2 V V e / / 



+ e/3 da;2 



+ e/3 9x2 



.2,00 _,.o , 5"2(^^^o)(^_^^) 



^ 1 + e/3 9x2 
(a;i,0) 



e2 92^0 



2 9X2 



7(1-X2) + -^(1-X2) 

1 + ep 



2 9^f(^-«) 



72^2 



We define the average wrt the fast variable in the horizontal direction: 

v{x) = — v{xi+y,X2)dy, e H\n'). 
27re Jo 



Then, one can see easly that for any x in fJ*^ 



1,00 , I r\\ ^^ \ 1 

1 + e/3 0x2 



2,00 I ^''^2 I n\ /I \ 

1 + e/? 0x2 



2 



^a;i,0) 



7(l-:r2) + -^(l-X2) 
1 + ep 



This means that the averaging process cancels the oscilations providing only macro- 
scopic terms still depending on e. Moreover one has the following compact form of 
the full boundary layer correctors 



0x1 



2Q2 2 (4-1) 



At this point, if one computes the boundary value problem that and 11} solve in 
the smooth domain, wc obtain the two following Robin and Wcntzel type problems. 
Namely, solves : 

- Aii^ = C, Va; e 17°, 

= V.T e r", = 0, Vx e r\ (4.2) 

(7X2 

is xi — periodic on Fin U Tout, 
whose explicit solution reads : 



u\x) = ^^(xl-^-^) 



(4.3) 



while the second order wall-law satisfies the folowing boundary value problem 

- Au^ = C, Va; e 17°, 



^2^e;3^ + -7r__, Va;er°, (4.4) 
6'a:;2 2 

2 



0, Va; G F^, is xi — periodic on Fjn U Fq 



4.2. Existence and uniqueness of the second order wall-law. Because 
problem (j4.4p contains second order normal derivatives as components of the boundary 
condition, (in the literature this kind of boundary conditions are called of Wentzell 
boundary conditions) the existence and uniqueness is not a standard result. Here we 
provide it. First we transform the second-order normal boundary term in a tangential 
term of the same order. Then using the appropriate test function space, we can apply 
Green's formula on tangential directions and symmetrise the bilinear form associated 
to the problem. 

Lemma 4.1. Under hypotheses \2. li the system (|4.4p admits a unique solution in 
if^'^(ri°) = {v e _ffpi(f7°); V e _ff^(F°)}, where H^i is the set of functions belonging 
to H^{ff), xi — periodic on T\a U Tout cmd vanishing on F^. 
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Proof. The boundary condition shall be transformed thanks to the first equation of 
(113) into 



2 dxn 



dx2 



Because P*^ does not intersect F, and thanks to the maximum principle, /3 > a.e. in 
Z+ U P. This implies that /3 > which allows the weak formulatiuon [TU] : 



1 7 
^ ^ ' 2/3 

7 



|i.|(x„0) 



-1 xi—L 



^ f du dv 



2-1=0 



2(3 \dxi dxi /po 



(C,w)oo - e-=(C,i;)ro, 
2(3 



where the third term of the Ihs vanishes thanks to the periodicity of the solution 
and of the corresponding test functions of H^{Vf'). We have obtained a symmetric 
problem. Because 7 € [— 1, 0[ and (5 g]0, 1], setting 



1 T / du dv 



e/3 



2/3 \dxi dxi / po 



l{v) = {C,v)n~e^{C,v)ro 
Zp 



one obtains a variational formulation where a is coercive, H^^{ilP) being endowed 
with the norm : 



l^llffi 



(pO). 



Moreover, a and I are continuous on H^^l^lP), thus the problem is solvable by the 
Lax-Milgram theorem. By the way, we derive the following energy estimates that 
describe the dependence of various norms upon e : 



I|m|Il2(po^ < V^c, 



du 



dxi 



< 



C 



L2(ro) 



Note that when e goes to zero, our approximation leaves i?^'^(f2'^) moving to 
iuro(^°)- loose the control over the tangential derivative on the boundary. 

□ 

In the particular case of a straight domain fi^ this unique solution is explicit and 
reads 



X2{l + e^^) e(/3-e7) 



l + e/3 



l + e/3 



(4.5) 



4.3. Macroscopic error estimate. When replacing the Poiseuille profile in fi" 
by u^ or u^, one can compute the corresponding error estimates. 
Proposition 5. Let he the solution of (|2.2p and u^ (resp. u^) be the solution of 
(|4.2|) (resp. (|4.4|) ). then Under hvvotheses [2A[ 



^u"^ — u 11^2(5^0-) < C'e^ , and \\u'^ — u ||^2(qo) < Ce^ . 
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Proof. Wc only compute the error of the second order approximation, the case of 
is identical. We take advantage of estimates obtained in proposition [2] by inserting 
the full boundary layer corrector between u*^ and : 



where we used the compact form exhibited in (|4.ip . Then, one gets 



Thanks to proposition [21 and the last estimate in the claim of theorem 13.11 one gets 
the desired result. □ 

Remark 4.1. This result is crucial: it shows that the oscillations of the first order 
boundary layer edu^/dx2{f3 — (3) are larger than the second order macroscopic contri- 
bution. It is also optimal (see section\^for a numerical evidence). This observation 
motivates the sections below. 

5. Multi-scale wall-laws. In this section we continue the investigation in the 
sense introduced above. We aim to compute a solution that exists in fl" as or 
but that performs a better approximation of the exact solution u"^ restriced to . 
Below we shall show that this concept provides some new multi-scale wall-laws. 

5.1. The first order explicit wall- law . How can first order correction be im- 
proved if the non-oscillating second order extension of Saffman-Joseph's condition 
does not help. The aswer below will be to take into account some multi-scale features. 
If we consider the full boundary layer corrector u\'°° , it solves p.Sp . Moreover, on 
the fictitious boundary , its value is easily computed, namely 



= e-^(xi,0)/3(xi,0). 



We use this value as a non-homogenous Dirichlet boundary condition on for a 
Poisson problem that is nevertheless homogeneous on F^. Indeed, we consider the 
following problem 

- AZY, = c, Vx e 

Z^. = e|^(x-i,0)/3(^,0), Vxer", (5.1) 

c on Fin U Tout, 
and we claim the following 

Proposition 6. Under hvvotheses \2.1[ one gets the following error estimates 



Proof. Following the same lines as in the proof of proposition [5l one inserts the full 
boundary layer approximation error r^'°^ u'^ ~ u\'^ : 

r^i ^u'- + ^U,:^ r^'°° - [U, - u\'°°] =: r^'°° - J. 
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The first part of the rhs has already been estimated (prop. [2]). It remains to estimate 
the last term J, that solves the following system : 



-AJ = 0, Vxen^, 
J = 0, Vx e r°, 

Xi 1 
e e 

J is Xi — periodic on Fin U F, 



J = eg(x„0)(/3 



-/3 



Vx e F^ 



Using a ?/2-linear lift s that takes away the F^ boundary term (which is exponentially 
small wrt e), and thanks to the Poincare inequality, we obtain 



l^ll 



< cial 



< cue 



where C13 and C14 are constants independent on e. □ 

Remark 5.1. The error in 0{e^) is only due to the first order boundary layer ap- 
proximation. Indeed the extension of the Poiseuille flow is only linear inside \ fi". 
Nevertheless, we avoid errors when neglecting microscopic oscillations in our macro- 
scopic problem as it was the case for and . 

5.2. A second order explicit wall-law. Extending the same ideas as in the 
subsection above, one sets the following multi-scale problem: find Ve G H^{Vl^) such 
that 



AVe = C, Vx e 17°, 



2 dxl 



(^,0), VxeF°, 



(5.2) 



^ Vc = 0, Vx e F^, Vc is Xi — periodic on Fjn U Fouti 
for which we can prove 

Proposition 7. Under hypotheses [KJl one gets 

II"' - ^£11^2(^0) < cise-^ 
where the constant C15 is independent on e. 

5.3. First order implicit wall-laws. Note that the standard averaged wall- 
laws u^,u^ are building blocks of explicit multi-scale approximations U^^V^ solving 
problems (|5.1|5.2p . In this part we look for an implicit approximation that avoids the 
computation of these lower order approximations. Indeed, at first order we propose 
to solve : 



-AT, 



c, Vx G 



e/?(^,0)^, VxGF°, 
e 0x2 

0, V.X G F^, Tj is xi — periodic on Fin U Fq 



(5.3) 



We give here a first result of this kind : 

Theorem 5.1. Under hypotheses \2.1l there exists a unique solution Te G H^i{Q°) 
of problem (j5.3|) . Moreover, one gets : 
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where Cig is a constant independent of e. 
Proof. There exists a unique solution Tg solving 
the weak formulation of (|5.3p reads : 



Indeed, under hypotheses 12.11 



a{u,v) := (Vu, ViOfjo + ( 



po 



(Vu, Vu)oo 



At the mieroscopic level, we suppose that does not eross F, thus there exists a 
minimal distance S > separating them. By the maximum principle, /3 is bounded: 
(3 £ [S; 1]. Thus 1/(3 is bounded a.e. The bilinear form a is continuous coercive in 
H^iiJlP), the linear form / is continuous as well, thus existence and uniqueness follow 
by the Lax-Milgram theorem. To estimate this new approximation's convergence rate 
we add and substraet Ue, the explicit wall-law between w*^ and T^. 



'bi 



(5.4) 



Q is the solution of the boundary value problem reading : 

-Ae = 0, Va; e O, 

dx2 0X2 

6 = 0, Va; G r^, is a;i — periodic on Fin U Tout- 



e = e/3 



V.T e r". 



We reexpress the boundary condition on F'' introducing a Robin like condition, 
namely : 



ae 

e - e(3— = e(3 

OX2 



au 
dx2 



8X2 



V.T e r°. 



(5.5) 



where the rhs is explicitly known. We have the following weak formulation : 
-{Ae,v)no ^-(^, v] +(Ve,Vi;)no =0, \fveH^^{n°), 



po 



where the space -ffpi(f7") contains i7^(il") functions vanishing on T^. Then using 
(15.51) one writes : 



e 



(e,tO = (ve,v.),,o+ = 



e(3' 



po 



dx2 0X2 ' 



po 



We remark that the rhs is in fact a boundary term of another comparison problem 



and we set z 



Lie where z is harmonic and solves 



=-(Az,t.)no-(Vz,Vw)no, e H^i{n" 



po 



Estimates of the gradient. We have recovered a simpler problem that reads 

a(e, v) = -(Vz, Vw)oo, Vu e H^^ (17°). 
Thanks to proposition [2] and proposition [6l one gets 



l!ve||^.(,,o) < ||Vz||^.(^o) < ||v(u'-^ 

where ii' is a constant independent of e. 



,0, < 2ci7e, 
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Estimate of the trace. The control on the interior term enables to recover trace 
estimates 

\mUr") < miL^^D I J^r^dx, < e||Ve||^.(,,„)||Vz||^.(,,o) < cj^e^ 

Final estimate. By the dual problem, and trace estimates above, we finaly obtain 

||6|Il2(oo) < ci8||e||^2(ro) < cwei, 
Recalling relation (|5.4p . one gets : 

1^1^1,11^2(00) < lkbi|lL2(no) + l|0|lL2(no)' 
which ends the proof. □ 

Remark 5.2. A similar implicit approach could be considered at second order. This 
should lead to consider a multi-scale Wentzel condition. It is an open problem to show 
existence, uniqueness and error estimates as in theorem \5.1\ in this case. 

6. Numerical evidence. We compute u^, a numerical approximation of the 
rough problem (|2.2p on the whole domain 17*^, e taking a given range of values in 
[0.1, 1]. Then, we restrict the computational domain to , and compute macroscopic 
approximations w^, u^,Z^£_a, Ve_A, T^^A, again for each value of e. We evaluate the 
errors w.r.t. Ua interpolating the latter exact solution over the meshes of the former 
ones. 

Computational setting. For every simulation, we use a P2 Lagrange finite element 
method implemented in the C++ code rheole43 [H]. Our computational domain is 
a channel of length L = 10 and of height h = 1. We assume a rough periodic bottom 
boundary defined by formula (|2.ip with 

(l + cos(yi)) 
7(2/1) := 2 ' 

where (5 is a positive constant set to 5e — 2. 

The rough solution u^. We compute over a single macroscopic cell x G w*^ := 
{xi € [0, 27re] and X2 S [/(xi/e),l]} and we assume periodic boundary conditions 
at {x2 = 0} U {xi = 27re}. For each fixed e, we mesh the domain lo'^ while keeping 
approximately the same number of vertices in the xi direction. This forces the mesh to 
get finer in the X2 direction in order to preserve the ratio between the inner and outer 
radius of each triangular element. With such a technique we avoid discretizations that 
could be of the same order as e. 

Cell problems. In order to extract fruitful information for macroscopic wall-laws, we 
compute first and second order cell problems. Again we impose yi-periodic boundary 
conditions. We truncate the upper infinite part of the domain by imposing a homo- 
geneous Neumann boundary condition at ?/2 = 10 after verifying that a variation of 
the domains height no more affects the results. In [T3], the authors show an expo- 
nential convergence w.r.t. to the height of the truncated upper domain towards the 
y2-infinite j/i-periodic cell problems (|3.2p and (|3.6p . this validates our approach. Cell 
problems are computed over a mesh containing (9211 elements and 4738 vertices). 



^http: //Ijk. Imag. f r/membres/Pierre ■ Saramito/rheolef / 
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figure 6.1: The traces /3a (2/1, 0) and7A(2/i, 0) 



eA 


Ua - "a 


"a - "a 


"a - "a 


Ua - ^e,A 


"a - ^e,A 


"a - T^£,A 


a 


l.ll 


1.4786 


1.3931 


1.768 


2-3.6 


1.6227 



Table 6.1 

Numerical orders of convergence for various approximations 



We extract solutions' trace on the fictitious interface T for both first and second order 
cell problems (cf. fig. 16. ip . and compute the averages /3 = 0.43215 and 7 = 0.29795. 
Macroscopic approximations: Classical & new wall-laws We compute the classical 
macroscopic wall- laws over uj^ ~ {x G lo'^ / X2 > 0}, a single periodicity cell of ^}P. 
We follow the same rate of refinement as described above. Then, we solve problems 
(14.214.41) . 

In the same spirit, we use both averages {(3,j) and oscillating functions 
/3(^, 0), 7(^, 0) as a non-homogenous Dirichlet boundary condition over the macro- 
scopic domain when solving (|5.ip and (j5.2p . To provide values at the boundary we 
use a Pi interpolation of the data extracted from the cell problems. 
For the implicit multi-scale wall-law, we solve system (j5.3p using the inverse of 
f^Aixi/e, 0) as a weight in the boundary integrals of the discrete variational formula- 
tion. 

Results. We plot fig. 16.21 the L^^ff') error computed respectively for approxima- 
tions presented above: — u^, — u^, — u^, — Z^e.A, W^a ~ ^c,Ai ""a ~ "^e.A- 
If we set CA = Ce", table [6Tl gives approximate numeric values of convergence rates. 
Interpretation. A first important result, visible fig. 16. 2i is that there is no difference 
between first and second order macroscopic wall-laws and u^. This proves that our 
estimates are actually optimal. It explains also why one could never distinguish first 
from second order approximations in [21 H] . 

Next, we remark that convergence orders are not better than those predicted by the 
estimates for u\^u\M^,a, while the error displayed for Ve.A is limited by the P2 
interpolation. Indeed, the H^{QP) error is of order 3 on the vertices but is worse 
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€ 



figure 6.2: L'^{n'^) error computed versus e 

elsewhere inside the elements. Nevertheless, the error — Vc^a is more than one 
order smaller than for uj^^U^^A for every fixed e. 

The fully explicit oscillating wall-laws Uc^A,^e,A provide better results than the im- 
plicit ones, u^jU^ and T,;. Indeed, in the former the shear rate du'^ /dx2{xi,0) and 
the second order derivative d^u'^ /dx2{xi,0) of the limit Poiseuille profile are explicit 
and included in the boundary condition, whereas the latter approximate this informa- 
tion as well. This leads to supplementary errors on the macroscopic scale for implicit 
wall-laws. 
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Appendix A. The cell problems. 

A.l. Various properties of the first order cell problem's solution. Ex- 
istence and uniqueness of solutions of system (|3.2p , liave been partially proven in [2] . 
The authors consider a truncated domain supplied with a non-local "transparency" 
condition, the latter is obtained via the fourier transform. We give here a rigorous 
proof in the unbounded domain framework. 

Proof, [of theorem 13. 1| In what follows we express the cell problem as an inverse 
Steklov-Poincare problem solved on the fictitious interface F. This allows us to char- 
acterize /3 the solution of (|3.2[1 on domains and P separately, as depending only 
on T], the trace on F. Wc apply domain decomposition techniques [20]. In a first 
step we give a simple proof of existence that guarantees the existence of the gradient 
in U P). The solutions of the cell problems arc not in the classical Sobolev 

spaces because the domain Z'^ is unbounded in the y2 direction: the solutions are 
only locally integrable. For this purpose, we define, for an arbitrary open set uj, 

D'^'Piuj) = {v G Ll^{uj)/D''v e hP{uj), Ma G = n, v is yi - periodic }. 

In the particular case when n ~ 1 and p = 2, we define Dq''^{uj) := {v £ 
D^''^{uj)/ v\gi^ = 0} := Vo(uj), which is a Hilbert space for the norm of the gradi- 
ent. 
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Lemma A.l. Problem (|3.2p admits a unique solution (3 belonging to D^'^{Z^ U P). 
Proof, [of lemma [XT] Wc define the lift s ~ 2/2 X[P] that belongs to U P). 

Setting jS ~ j3 — s. the lifted problem becomes 

'-5t, inZ+UP, 
on P°, 

where (5r is the dirac mesure that eoncentrates on the fictitious interface F. The 
equivalent variational form of this problem reads 

a{P, v) = l{v), V G dI-^{Z+ U P), (A.l) 

where a{u,v) = (Vu, Vu)2+up ^iiid l(v) = — (Vs,Vu)p. These forms are obviously 
continuous bilinear (resp. linear) on i:»o'^(Z+UP)xi:)o'^(Z+UP) (resp. dI-^{Z+UP)). 
Because of the homogeneous boundary condition the semi-norm of the gradient is a 
norm. By Lax-Milgram theorem, the desired result follows. □ 
We define the following spaces 

Vi = D1'2(Z+), V2 = {v G H^{P) s.t. v\po =0, V isyi- periodic } 
Vi,Q = {veVi, v\r = 0}, V2,o = {veV2, v\r ^ 0} 

{ri e s.t. r/ = v\r for a suitable v G dI'^{Z+ U P)}. 

Lemma A. 2. The following domain decomposition problem is equivalent to (jA.ip ; 
we look for /?2) G Vi x V2 such that 



ai(/3i,w) (V/3i,V«)z+ =0, V« G V^i,o, 
/3i = /32, OK F, 

a2{f32,v) (V/3i, Vw)p = -(Vs,Vw)p = 0, Vw G ^2,0, 

a2(/?2,7^2/i) = -(Vs,V7^2A^) -al(/3l,7^lAi), V^i G A, 



(A.2) 



where TZi denotes any possible extension operator from F to Vi. 

Proof, [of lemma IA.2| Let us start by considering the solution /? of (|A.1[) . Setting 
f3i = P\z+,l32 = /3|p, we have that (3, G Vi and that (|X2|) .l. (|X2l) .2 and (|X2)) .3 
are trivially satisfied. Moreover, for each /i G A, the function TZji defined as TZ^ = 
T^itiXz+ + T^2lJ'XP belongs to Vq. Therefore we have a{P,TZfi) ~ {f,TZfi)yfi G A 
which is equivalent to (|A.2[) .4. 

On the other hand, let Pi be the solution of (jA.2p . Setting /3 — PiX[z+] + P2XIP] from 
(|X2|) .2. it follows that V/3 G L'^{Z+ U P), and P\po = 0. Then taking v € Vq we set 
^ = i'|r G A. Define TZfi as before; clearly (u,; — TZi^L) G V^^o and from (|A.2p .l. (|A.2p .3. 
(|X2|) .4 it follows that 

a{P,v) = ^[ai(/3j,Ui - 7^^^) + aj(/3i, 7?.i/i)] = -(V.s, V7?.2Ai)p 

■i 

= -(l,//)r = -{1.v)t = -(Vs, V?;)p. 

□ 

T/ie Steklov-Poincare operator. The Steklov-Poincarc operator S acts between the 
space of trace functions A and its dual. More precisely, applying Green's formula and 
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setting HiT] to be the harmonic lift in Z+ (rcsp. P) for all ?y G A, we have 



where < •, • > denotes the duality pairing between A' and A. In particular, taking 
TZifi ~ HilJ-, we obtain the following variational representation : 



< Sr],fi>=^^ai{Hir],Hifi), V77, ^ G A. 



The linear form on A. We set l(fi) as follows : 

= -(VS, WH2fi)p = (1, -^H2li ] = (1, ^i)T■ 



, dx2 J p 
Lemma A. 3. The problem: 

find ?/ e A such that < Sjj,j-i >= ^fJ. G A, (A. 3) 

admits a unique solution. Moreover this is equivalent to solve (jA.2[) . 

Proof, [of lemma IA.3| We use the Lax-Milgram framework : 

- Continuity : 

< STJ,fI>< |lVi77]|1^2(2+up)|lVffA«llL2(Z+UP) ^ C2ohllAllMllA, 

by well know estimates for solutions of elliptic boundary value problems |15j . 
For H2 this can be computed explicitly (see below). The continuity of I is 
obvious. 

- Coercivity 

< Stj,7j >= \\VHr]\\l2^z+uP) ^ C2i||i/2f7||^i(p) > C22||?7lli- 

Then applying Lax-Milgram theorem one gets the desired result. 
To prove the equivalence between (|A.3p and (|A.2p . it suffices to separate the har- 
monic lift Hi and the solutions of the Poisson problem with homogeneous boundary 
conditions and the result follows as in [20] p. 10. □ 

The harmonic extension in Z'^ named Hi We set ry € A. By decomposing in j/i-fourier 
modes, one gets that the solution of : 

|A^ = 0, VyeZ*, ^^^^ 

[ /3 = 77, VyeT, 

rewritten as /? = J2k I^k{y2)e^''^^ , Vy G Z+ should satisfy the following system of 
ODE'S : 

/3^' - fc2/3fc = 0, 2/2 e M+ 
Pk{0) = Vk, 2/2 = 
/3fc(2/2)GL°°(R+;C), 
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where rjk = Jq^ e ^^^^ '>l[yi)dyi are 7?'s fourier coefBcients on F. The solution Pz+ is 
expUcit and reads 

oo 

^^i'? = /?|z+ = E VyGZ+. (A.5) 

k—~oo 

To show exponential convergence towards zero of /3 — /? and V/? when 2/2 — > 0, we use 
the same arguments as in the second part of [3], theorem 2.2.1 p. 637, whose proof is 
omitted. 

Proposition 8. There exists a\ > (47r)^/9 such that the solution of problem p.2p 
satisfies 

~ ^llL3(Z+UP,e°i«2) ^ C23||V/3||^2(2+Ljp^eai„2) < C24, 

which implies also (3's and V/3 's exponential decay in the 7/2 direction. 
□ 

A. 2. The second order boundary layer. Proof, [of proposition [3] Problem 
p.6p is equivalent to solve : 

rA7-2x[P], VyeZ+UP, 
[7 = 0, Vy G P°. 

This, under the previous domain decomposition form, reads: find (7z+,7p) such that 

' (V72+,Vz;)z+ =0, V«eifi(Z+), 
lz+ = IP, on r, 

(V7P,Vz;)p - -(2,w)p, yveH^^poiP), ^ 
^ (V7P, V7^p/i)p = -(2,7^PAi)p - (V7^+,7^z+/^)z+, V^i e H^T). 

Following the same lines as the proof above, we write the interface problem : 
< SX,^i>={VHpX,VHpfi) + {VHz+X,VHz+fJ.), V^i e H^(T), 

= -(2,i7pAi) - {vg■2,VHp^_i) =: i{ji), y^i e h^t), 

where Q2 is the solution of the homogeneous Poisson problem : 

A^a = 2, yye P, 
Qi = 0, Vy G P" U F, 
Qi is yi — periodic. 

One gets the continuity of the linear form again, thanks to the properties of the 
harmonic lifts [121 IH] '■ 

\l{ll)\ - |-(2, i/p/i) - (Vg2,Vi/pAi)| < C25|li/p/i|lHl(p) < C26|1a^|1^1(p)- 

And again, by the Lax-Milgram theorem, one gets the desired result. □ 
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